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We consider scattering atoms from a fully Bose-Einstein condensed gas. If we take these atoms 
to be identical to those in the Bose-Einstein condensate, this scattering process is to a large extent 
analogous to Andreev reflection from the interface between a superconducting and a normal metal. 
We determine the scattering wave function both in the absence and the presence of a vortex. Our 
results show a qualitative difference between these two cases that can be understood as due to an 
Aharonov-Bohm effect. It leads to the possibility to experimentally detect and study vortices in this 
way. 



The achievement of Bose-Einstein condensation in 
trapped alkali vapors has generated renewed theoretical 
and experimental interest in the subject of quantized vor- 
tices, since the quantization of circulation, as well as the 
presence of persistent currents, are amongst the most 
striking features of superfluidity in a Bose-condensed 
fluid. In the case of superfluid helium, vortices have been 
studied extensively 0-^] . Because of the strongly inter- 
acting nature of such a liquid however, the understand- 
ing of these phenomena is mainly phenomenological. In 
contrast, in weakly- interacting Bose condensed gases the 
theoretical study of vortices can be based on a micro- 
scopic theory. 

Starting from such a microscopic viewpoint, problems 
that have been studied in detail are the stability of a vor- 
tex in a harmonic trapping potential, the related ques- 
tion of its experimental creation, and the detection of 
these vortices |5Hl9fl. The creation of a stable vortex 
appears to be possible by rotating an external trapping 
potential that is sufficiently anisotropic with a frequency 
larger than some critical value In this case, a vortex 
configuration becomes the thermodynamic equilibrium of 
the system. This equilibrium can be reached in two ways, 
either by first cooling the gas through the transition tem- 
perature and then rotating the trapping potential, or by 
first rotating and then cooling the gas. The latter seems 
to be the most favorable experimentally ||. Concerning 
the detection of a vortex, several ideas have been put 
forward. It might be possible to detect the core of the 
vortex after the Bose-condensed cloud is released from 
the trap and has been allowed to expand for some time 
|pT| . Also, the interference fringes between two conden- 
sates containing vortices will exhibit dislocations Jl2],[l3| . 
In addition, the excitation frequencies of the collective 
modes of the condensate are shifted due to the presence 
of a vortex (TJdT^] . Finally, it has been proposed to study 
vortices by optical means Jl7| ]. 

In addition to these ideas, in this paper we suggest and 
investigate a method which is based on the scattering 
of identical atoms off the Bose-condensed cloud, analo- 



gous to the occurrence of Andreev reflection |2(j] from a 
superconductor-normal metal interface J2l|] . In the latter 
case, the scattered electrons probe both the magnitude 
and the phase of the superconducting order parameter. A 
similar experiment with a trapped Bose condensate will, 
besides the condensate density, also probe gradients in 
the phase of the Bose-condensed cloud. This is desirable, 
since the vortex core is typically small compared to the 
overall size of the condensate, whereas the superfluid ve- 
locity due to the presence of the vortex extends over the 
entire Bose condensed cloud. Because the velocity field of 
a condensate with and without a vortex is qualitatively 
different, the scattering of atoms might be a possible way 
to detect a vortex. 

Note that a similar idea has been put forward recently 
in the context of superfluid 4 He, to study the quantum 
sticking, scattering and transmission of 4 He atoms from 
superfluid 4 He slabs |^2[ . In an experiment with a Bose- 
condensate, the energy of the incoming atoms has to be 
rather small to see an effect, namely of the order of the 
chemical potential of the condensate. A convenient way 
to realize this experimentally is by using optical means to 
extract the atoms that are to be scattered from the con- 
densate itself. After displacing these atoms with respect 
to the condensed cloud, they are subsequently released 
and measured by the usual time-of-flight measurements 
or in situ optical imaging [^3) . 

Bearing this in mind, we consider the whole scattering 
process to take place in an external confining potential. 
Such a process can be described by the time-dependent 
Bogoliubov-deGennes equations. To lowest order approx- 
imation, the condensate density profile can be taken to 
be static, which implies that the time dependence of the 
mean-field interaction due to the possible excitation of 
one or more collective modes is neglected. In the static 
approximation, the coherence factors /(r,i) and g(r,i) 
for the atom satisfy 
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Here, T 2S = Airfi 2 a/m is the s-wave approximation 
for the two-body scattering matrix and V ext (r) — 
^.mwfx 2 /2 denotes the external trapping potential. 
The inter-atomic scattering length is given by a. Fur- 
thermore, the static density profile of the condensate is 
denoted by p(r) and its phase by x( r )- Outside the con- 
densate, g(r, t) has to be zero and Eq. (|l|) reduces to the 
Schrodinger equation for an atom in a harmonic poten- 
tial. 

Of course, the colliding atom can in principle excite the 
condensate. To include the possibility of exciting collec- 
tive modes, we simply have to replace p(r) and x( r ) by 
their operator counterparts, p(r,t) and x(r,t), respec- 
tively. The latter two can be expanded in terms of the 
creation and annihilation operators b\ and hi of the col- 
lective modes with mode functions Ui(r, t) and Wj(r, t) as 



p(r,t)=p{r) + y/p{r) ^ [«i ( r , *) ~ «i 0, *)] 
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that |0) can in principle be an arbitrary solution of the 
Gross-Pitaevskii equation and can in particular also be a 
vortex or a kink solution. 

In the remainder of this work, we apply the static 
approximation and include only the ground state |0). 
The resulting equations describing our scattering prob- 
lem then reduce to Eq. (Q) . This is expected to be accu- 
rate in the case of the scattering of a single atom at ener- 
gies several times larger than the chemical potential. At 
these energies, the kinetic energy of an incoming particle 
is still relatively large and the overlap of its wavefunction 
with that of the low-lying collective modes small. In ad- 
dition, we take the external potential to be rotationally 
symmetric around the z-axis, i.e., to x — lu v = uj r , and ap- 
proximate the cigar-shaped traps used in experiments by 
assuming translational invariance in the z-direction. It is 
then convenient to write the Laplace operator in cylindri- 
cal coordinates (r, z, 0), i.e., V 2 = 
and expand /(r,t) and g(r,t) as 
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and 
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Inserting this expansion into Eq. ([!]), the Bogoliubov- 
deGennes equations for the wavefunctions f n (r, t) and 
(3b) g n (r, t) become, 
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If we denote by \a) the number states in the Fock-space 
spanned by the creation and annihilation operators, the 
different channels of our scattering problem are given by Here; the matrix elements H nn > and V nn > equal 
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and the equations of motion for these different channels 
read 
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Again, outside the condensate g a (r,t) has to be zero. It 
is clear that if the condensate is initially in the ground 
state |0), after colliding with the incoming atom it can 
be in the ground state or any of the excited states. Note 
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Note that we consider the solution of Eq. (Jl]) in the pres- 
ence of a vortex with a winding number I which is aligned 
in the direction of the z-axis. This implies that the con- 
densate wavefunction is given by ip(r, t) = p(r, z,t)e l1 ^ . 

Upon solving the resulting equations we treat the con- 
densate in the Thomas-Fermi approximation, both in the 
presence and in the absence of a vortex. Thus, the den- 
sity profile and the phase are given by 
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and 
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where O [x] denotes the Heavy-side step-function. To 
lowest order, the change in the chemical potential due 
to the presence of a vortex can be neglected Jl6| . 

The method we use to calculate numerically the time- 
evolution of an initial wave-packet describing the incom- 
ing atom is based on Cayley's finite difference represen- 
tation of the Schrodinger equation p3 , 



a iHAt 1 " At/2 
l + iHAt/2 
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Here, At is the time step in the discretized Schrodinger 
equation. This algorithm is second-order accurate in time 
and unitary, i.e., the norm of the wavefunction is con- 
served up to computer accuracy. 

The results of our numerical calculations are show in 
Fig. H and Fig. |[ corresponding to a condensate with- 
out and with a vortex, respectively. In these calcula- 
tions, the initial wavefunction of the incoming atom is 
given by the product of a Gaussian in both the r and 
the <f> direction, centered around r — tq and (/) — ?:. In 
principle, we could have plotted the probability distribu- 
tion P(r,cj),t) = |/(r, (j), t)| 2 — \g(r,cj),t)\ 2 as a function 
of r and <f>. As time evolves, the hole part g(r, <fi, t) of 
the atomic wavefunction becomes nonzero if the particle 
part f(r, (f>, t) starts to overlap with the condensate. Also, 
the probability distribution P(r, </>, t) develops nodes in 
the radial direction. The wavelength of these nodes is a 
measure of the kinetic energy of the incoming particle. 
However, to present our data more comprehensively, we 
integrate with respect to r and denote the resulting angu- 
lar probability distribution by P(4>, t). Note that P(<p,t) 
is normalized to one, and that this normalization is con- 
served numerically. The integrated distribution is shown 
at four different times, corresponding to t = 0, approxi- 
mately half a period of oscillation t = 3.5 , and two 
snapshots at t = 1.8 ui^ 1 and t = 2.0 u)^ 1 . 

There are three main differences between Fig. [l] and 
Fig. H First, in Fig. [I] there is scattering in the forward 
direction, whereas in Fig. || there is almost no probability 
for scattering in this direction. Second, in Fig. ^, there is 
almost no scattering in the backward direction, whereas 
in Fig. ^ there is scattering in this direction. Third, and 
most important, in Fig. ^ the scattering is symmetric 
around <ft — 7r, whereas in Fig. this symmetry is clearly 
broken. The first and second point are related to the 
presence of the vortex core, which causes the reflection 
in the forward direction to be decreased in favor of the 
backward direction. This has been checked numerically 
by inserting the core, but not the phase of the vortex into 
the Bogoliubov-deGennes equations and solving for the 
wavefunction of the scattering particle. The third point 
however, cannot be related to the presence of the core, 



which is also rotationally symmetric. Instead, it is related 
to the phase of the condensate, which in the presence of a 
vortex indeed breaks rotational invariance. To make this 
point more clear, we present here a qualitative argument 
based on a semi-classical calculation of the Berry-phase 
associated with an atom moving around the vortex. The 
calculation shows that there is an Aharonov-Bohm effect 
p5[ in both the particle and the hole parts of the wave- 
function, which causes the total atomic wavefunction to 
pick up a non-trivial phase factor. 



2.0 



T3 1.5 
CO 



o 



c 

13 



1.0 



~ 0.5 
Q_ 

0.0 





in] 




\ [3] / 


1 1 


/ 




\\ [2] / 


_ / ■ 

r\ 1 
\ / 



0.0 2.0 4.0 

§ (in units of Rad) 



6.0 



FIG. 1. The angular probability distribution P(<f>,i) at 
four different times, in the absence of a vortex. At t = 0, 
P{4>,t) is a Gaussian centered around 4> — ""j line [1]- Line [4] 
is P((f),t) after approximately half an oscillation, t — 3.5 ui^ 1 . 
The short dashed [2] line and the long dashed line [3] are snap- 
shots at intermediate times, t = 1.8 w" 1 and t — 2.0 c^T 1 , 
respectively. The calculation is performed for 85 Rb in a trap 
with u r /2% = 200#z. 
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FIG. 2. The angular probability distribution P{<f>, i) at 
four different times, in the presence of a vortex with 1 = 1. 
The labeling of the curves, and the times to which they cor- 
respond are the same as in Fig. m. 
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The Berry-phase 6 associated with the spin-like na- 
ture of a two-component wavefunction is the phase picked 
up when the atom moves adiabatically along a classical 
trajectory C p3). If the two-component wavefunction 
IV'(x)) describes the 'spin' degrees of freedom, it can be 
expressed as 

e = i f dx-(V>(x)|V|^(x)) . (12) 
Jc 

In the present case of the Bogoliubov-deGennes equa- 
tions, we thus get 

8 = 1 ( dx- [/* (x)V/(x) - ff*(x)V.g(x)] 
Jc 

Til I 

= -- / rfx.v s (x)[|/(x)| 2 + | 5 (x)| 2 ] . (13) 

Here, the superfluid velocity v s = Wx/2ro, and the co- 
herence factors /(x) and g(x) correspond to one of the 
local eigenfunctions of the Bogoliubov-deGennes equa- 
tions, i.e., without the kinetic term ?i 2 V 2 /2m. In the 
presence of a vortex located at the origin with winding 
number I, the Berry-phase for the special case of a circu- 
lar trajectory at radius r therefore becomes 

= -27rf[|/(r)| 2 + | fl (r)| 3 ] . (14) 

This indeed shows that there is a nontrivial phase accu- 
mulated when the atom moves around the vortex, which 
causes the interference pattern to become asymmetric 
with respect to cf> = tt. Note that it is essential that the 
particle and the hole components of the atomic wave- 
function contribute with a different sign to the integral 
in Eq. (fL3|). To understand the direction of the shift in 
the interference pattern, consider two interfering paths 
that enclose the vortex. Seen from <fr = tt, one is passing 
on the left side, and the other is passing on the right side 
of the vortex. Because the accumulated phase for the 
contour C = Cr — Cl is negative in the presence of a 
vortex with I = 1, we again get constructive interference 
if Cr, becomes shorter and Cr longer. Therefore the in- 
terference pattern for a particle coming from <j) — tt shifts 
towards </> = 0, as seen numerically. 

In summary, we have presented a possible method to 
detect the presence of a vortex. This is realized by collid- 
ing atoms with the Bose-condensed cloud that are iden- 
tical to the atoms in the condensate. The most impor- 
tant signature of the vortex is an angular asymmetry in 
the scattered wavefunction, which we have related to the 
Aharonov-Bohm phases picked up by the particle and 
hole parts of the wavefunction. We note that to use the 
detection mechanism in the way presented here, the trap 
should not rotate. Thus after creating a vortex, the rota- 
tion of the trap must be stopped, in which case the vortex 
in principle becomes unstable. However, recent studies 
have shown that the time needed for the vortex to dis- 
appear from the condensate diverges with the number of 



particles in the condensate at least as N^ 5 |l8| , |2"7| ] . Thus, 
for large No this presents no difficulty for our detection 
scheme. Moreover, a detailed picture of the dynamics of 
the unstable vortex might be obtain by means of multi- 
ple collisions. In this paper however, we considered only 
a time-independent vortex configuration, and the effect 
of vortex dynamics on the scattering wavefunction is a 
topic of future investigation. 
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